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Abstract —Since the seminal work of Grossglauser and Tse (B, 
the two-hop relay algorithm and its variants have heen attractive 
for mobile ad hoc networks (MANETs) due to their simplicity 
and efficiency. However, most literature assumed an infinite 
buffer size for each node, which is obviously not applicable 
to a realistic MANET. In this paper, we focus on the exact 
throughput capacity study of two-hop relay MANETs under the 
practical finite relay buffer scenario. The arrival process and 
departure process of the relay queue are fully characterized, 
and an ergodic Markov chain-based framework is also provided. 
With this framework, we obtain the limiting distribution of the 
relay queue and derive the throughput capacity under any relay 
buffer size. Extensive simulation results are provided to validate 
our theoretical framework and explore the relationship among 
the throughput capacity, the relay buffer size and the number of 
nodes. 

Index Terms —Mobile ad hoc networks, finite buffers, two-hop 
relay, throughput capacity, queuing analysis. 

I. Introduction 

The self-autonomous and inherent highly dynamic char¬ 
acteristics of mobile ad hoc networks (MANETs) make the 
routing schemes based on opportunistic transmission widely 
applied. Among them, the two-hop relay algorithm which is 
first proposed by Grossglauser and Tse JT) and its variants m, 
have become a class of attracting routing protocols due to their 
simplicity and efficiency 0. The basic idea of this protocol 
is that the source node can transmit packets to the destination 
directly, or to a node (served as a relay) it encounters, later 
the relay node carrying the packet forward it to the destination 
node. Thus, each packet travels at most two hops to reach the 
destination. 

By now, a lot of work has been committed to studying the 
throughput capacity of two-hop relay algorithm in MANETs 
IfTl , Il2l , PI - 1^ . Grossglauser and Tse (2002) |[T] studied 
the order sense of per node capacity and showed that it is 
possible to achieve a 0(1) throughput by employing the node 
mobility, which is a substantial improvement for the static 
networks considered by Gupta and Kumar (2000) JT). Neely 
et al. (2005) E) computed the exact capacity and the end- 
to-end queuing delay for a cell-partitioned MANET, and a 
fundamental tradeoff between throughput capacity and delay 
has been developed. Sharma et al. (2006) a proposed a 
global perspective for delay and capacity tradeoff in MANETs. 


They considered a general mobility model and related the 
nature of delay-capacity tradeoff to the nature of node motion. 
The throughput capacity with packet redundancy has been 
researched in 0, while the throughput capacity under power 
control has been examined in 0. For a detailed survey, please 
see 0 and the reference therein. 

It is notable that all available work mentioned above as¬ 
sumed the node buffer is infinite. However, this assumption 
will not hold for a realistic MANET obviously. In 0, J. 
Herdtner and E. Chong have studied the throughput and 
storage tradeoff, and they showed that the limited relay buffer 
will degrade the throughput capacity. Since they only provided 
a scaling law relationship, the exact throughput capacity of 
two-hop relay MANETs under finite relay buffer size remains 
largely unknown by now. 

As a first step towards this end, in this paper, we analytically 
study the exact throughput capacity of two-hop relay MANETs 
with the consideration that the relay buffer of each node, 
which is used for storing other nodes’ packets, is strictly 
bounded. The main contributions of this paper are summarized 
as follows. 

• Considering the soure-to-relay transmission under finite 
relay buffer scenario, we carefully compute the arrival 
rate at the relay queue. By utilizing the occupancy prob¬ 
ability technique, we exactly characterize the departure 
process of the relay queue. 

• Based on the queuing process of the relay node, a finite- 
state ergodic Markov chain model is constructed to obtain 
the limiting distribution of the relay queue. With this 
framework, we proceed to derive the exact throughput 
capacity under any relay buffer size. 

« Extensive simulation results are provided to validate our 
new analysis model and explore how the throughput 
capacity varies with the relay buffer size and the number 
of nodes. The results indicate that the throughput capacity 
under finite relay buffer cannot stay constant as the 
network size growing, which is quite different from the 
inhnite buffer scenario. 

The remainder of this paper is outlined as follows. The 
system models and a modihed two-hop relay scheme are 
introduced in Section m In Section mil we develop the 
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relay node and try to send a packet, if the relay queue is full, 
this transmission fails leading to the packet loss and the energy 
waste. To avoid this phenomenon, we introduce a handshake 
mechanism before each source-to-relay transmission to con- 
hrm the relay buffer occupancy state, if the relay queue is 
not full, the transmission can be conducted; else, the source 
node remains idle. At any time slot for each cell containing at 
least two nodes, the cell executes the modified two-hop relay 
(M2HR) algorithm which is summarized in Algorithm [T] 


Fig. 1. The local queue and relay queue in a node. 

ergodic Markov chain-based framework to fully characterize 
the queuing process of relay nodes and proceed to obtain the 
exact throughput capacity. The simulation results are provided 
in Section |Iy] Finally, we conclude this paper in Section [V] 

II. Preliminaries 

A. System Models 

Network model. The cell partitioned network model ID is 
adopted. The network is partitioned into C non-overlapping 
cells of equal size. N mobile nodes roam from cell to cell 
over the network according to the independent and identically 
distributed (i.i.d) mobility model m, m. With i.i.d mobility 
model, at the beginning of each time slot, each node selects 
a cell among the network uniformly and independently, then 
stays in this cell during this time slot. Thus, i.i.d mobility 
model can be served as the limit of infinite mobility. Further, 
we assume that time is slotted with a fixed length, and during 
each time slot, only one node in each cell can transfer exactly 
one packet to another node in the same cell. Nodes located 
in different cells cannot communicate with each other. We 
assume each node has a local queue and a relay queue. The 
local queue is used to store the self-generated packets and 
there is no constraint on it; while the relay queue is used to 
store the packets from other nodes and the buffer size is set to 
be B (packets) ||8l. The reason for this buffer assumption is 
elaborated in a, where the ingress buffer and the internal 
buffer correspond to the local buffer and the relay buffer 
respectively. Thus, the node model can be represented in 
Fig.IU 

Traffic model: We consider the traffic pattern widely adopted 
in previous studies 0, Eol, where N is even and the source- 
destination pairs are composed as follows: 1 -(A 2, 3 -fA 4, 
• • •, {N — 1) -H- N. Thus, there are in total N distinct 
unicast traffic flows, each node is the source of a traffic flow 
and meanwhile the destination of another traffic flow. The 
exogenous packet arrival at each node is a Bernoulli process 
with rate A packets/slot. 

B. A Modified Two-Hop Relay Algorithm 

In this section, we make a modification on the traditional 
two-hop relay algorithm, to be applicable under the hnite relay 
buffer scenario. We consider that a source node encounters a 


Algorithm 1 M2HR routing algorithm 
1: if There exist source-destination pairs within the cell then 
2: With equal probability, randomly choose such a pair to 

do source-to-destination transmission. 

3 : if The source contains a new packet then 

4 : The source transmits the new packet to the destina¬ 

tion. 

5 : else 

6: The source remains idle. 

7 : end if 

8 : else 

9 : With equal probability, randomly designate a node 

within the cell as the sender. 

10: Independently choose another node among the remain¬ 

ing nodes within the cell as the receiver. 

II: Flips an unbiased coin; 

12: if It is the head then 

13 : The sender conducts source-to-relay transmission 

with the receiver. 

14 : if The sender has a new packet, meanwhile the 

receiver informs that the relay queue is not full then 
15 : The sender transmits that packet to the receiver. 

16: else 

17 : The sender remains idle. 

18: end if 

19 : else 

20: The sender conducts relay-to-destination transmis¬ 

sion with the receiver. 

21: if The sender has packets destined for the receiver 

then 

22: The sender transmits a packet to the receiver. 

23 : else 

24 : The sender remains idle. 

25 : end if 

26 : end if 

27 : end if 


III. Throughput Capacity Analysis 
In this section, we hrst introduce some basic probabilities. 
Then, an irreducible ergodic Markov chain-based framework 
is established to explore the occupancy distribution on the 
relay queue. In order to solve the Markov chain, we further 
investigate the departure process of the relay queue. Based 
on this framework, we proceed to derive the exact throughput 
capacity. 





































A. Some Basic Probabilities 


For a given time slot and a particular cell, we denote by p 
and q the probability that there are at least two nodes and at 
least one source-destination pair in a cell, respectively. Then, 
the same as Q, we have 



Under the M2HR algorithm, we denote by psd. Par and 
Prd the probabilities that a given node has a chance to 
conduct a source-to-destination transmission, source-to-relay 
transmission, and relay-to-destination transmission at a given 
time slot, respectively. Then we have 


C 

Pad = 

Psr — Prd — 


C{p - q) 

2N 


(3) 

(4) 


The derivations of probabilities p^d. Par and prd are omitted 
here due to space limit, and please kindly refer to Appendix 
B in H for details. 


B. Analysis of the Local Queue and the Relay Queue 

Under the M2HR algorithm, each packet experiences at 
most two queuing processes, i.e., the packet dispatching pro¬ 
cess at the local queue and the packet forwarding process at 
the relay queue (if the packet cannot be transmitted to the 
destination directly). 

Local Queue: Due to the i.i.d mobility, the local queue can 
be represented as a Bernoulli/Bernoulli queue, where in every 
time slot a new packet will arrive with probability A, and a 
service opportunity will arise with a corresponding probability 
PS which is determined as 

Fs(A) = Pad + Par (1 “ Pb) , (5) 


where Pb denotes the probability that the relay queue is full. 
Note that the Bernoulli/Bernoulli queue is reversible im, so 
the output process is also a Bernoulli flow with rate A. 

Relay Queue: Since a specific packet from the output 
process of a local queue is transmitted to relay nodes with 
probability each of the N — 2 relay nodes are 

equally likely to encounter the source, and for each relay node 
there are in total N — 2 independent output processes of local 
queues may arrive at its relay queue. We denote by A the 
packet arrival rate at a relay queue when it is not full (when 
a relay queue is full, its input rate is 0), then we have 

A • (1 - Pb) + 0 • Pb = (iV - 2)A ■ - 2), 

Ms(A) 


=> A = 


Xpar 

PsW 


( 6 ) 


We denote by pR{k) that the service rate when the relay 
queue contains k packets and it is clear that pu{k) > 0 
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Fig. 2. State transition diagram for the relay queue occupancy process. 


for 1 < k < B. Thus, the relay queue can be modeled 
as a discrete Markov chain as illustrated in Fig. |2] and the 
corresponding one-step transition matrix of the relay queue 
occupancy process is given by 


1-A A 


P = 


PR{k) l-X-pR{k) X 


Pr{B) 1-pr{B) _ 

(7) 

For this Markov chain we have the following lemma. 

Lemma 1: The limiting distribution of the relay queue 
occupancy process exists and is unique, and is equal to the 
stationary distribution. 

Proof: For any state i,j € S , i < j , there exist m and 
m, such that p[p > 0 and > 0. For example, p [ j ~^'^ = 
Xb-i) > 0 andp^7-*) = pR^j) x pR^j-1) x ■ ■ ■ x pR{i + l) > 
0. Thus, any two states i,jGS can communicate with each 
other, which is denoted by i •<->• j . This type of Markov chain 
is called irreducible (see, for example. Definition 4.1 in IfT^ 
or [13, ch. 4, p. 168]), and in this Markov chain all the states 
are in the same class. From 113 it follows that all the states 
are recurrent. 

We denote by d{i) that the period of state i. Note that pgg^ = 
1-A, p^Jl = 1 - pr{B), and = 1 - A - pR{k), for 
0 < k < B. Thus, for any state i G S, we have d{i) = 1. 
The state i and the Markov chain are called aperiodic 113 . 
Since all the states are recurrent and aperiodic, the relay queue 
occupancy process is an irreducible ergodic Markov chain. 
Referring to [12, ch. 5], the limiting distribution of the relay 
queue occupancy process exists and is unique, and is equal to 
the stationary distribution. ■ 

We use n = {7r(0), 7r(l), • • •7r(P)} to denote the limiting 
distribution of the relay queue. By lemma [T] we have 


n ■ P = n, 


A7r(0) = /iB(l)7r(l), 
A7r(l) = pR{2)Tr{2), 


A7r(P - 1) = pr{B)tt{B). 


(8) 
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Fig. 3. Examples of occupancy. 


Combining with the normalization equation 

the limiting distribution of the relay queue is given by 

-1 

\k \ 

r(0) = 


r{k) = 



k'Riky. 


PR{k)l) 

(9) 


(10) 

PR{k) X pR{k-l) X • 



From the ® and (fTOl i. we can see that in order to derive the 
limiting distribution of the relay queue, we need to compute 
the service rate pij(fc) when the relay queue is on state k. 


C. Computation of the Service Rate at Relay Queue 

il) 

We denote by p), ’ the probability that the relay queue 
has k packets and these packets are destined for i different 
destination nodes, 1 < z < fc. Due to the i.i.d mobility model, 
we have 


From the O, it is clear that in order to compute pij(fc), 
we need to derive To address this issue, we utilize 
the occupancy technique (see, for example, [14, ch. 1]). We 
represent the packets by ’stars’ and the N — 2 destination 
nodes by the spaces between N — 1 ’bars’. For example. Fig. [3 
represents one packet is destined for the first node, no packet 
is destined for the second and third nodes and two packets are 
destined for the fourth. 

Notice that N — 2 nodes need N — 1 ’bars’, since the first 
and last symbols must be ’bars’, only TV — 3 ’bars’ and k 
’stars’ can appear in any order. The number of all possible 
permutations is given by 


Ek 


i^-n 


( 12 ) 


Considering these packets are destined for i different destina- 

(i) 

tion nodes. The number of permutations E^, is given by 
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i — 1 + k — i 
k — i 


According to the classical probability we have 
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Substituting ([l4li into (ITTIl we have 
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Vi/' \k-iJ 

2^—pra+y- 


. ) 


TV- 2 


(13) 


(14) 
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D. Throughput Capacity 

Definition: Throughput Capacity: For a MANET under 
M2HR and with a given packet arrival rate A, the network is 
called stable if the queue length in each node (thus the average 
delay) is bounded. The throughput capacity Tc of the network 
is then defined as the maximum value of A the network can 
stably support. 

Theorem 1: For a MANET with finite relay buffers, the 
throughput capacity under M2HR satisfies 

Tc = max{X : A < p,s{X)}- (16) 

Proof: It is notable that Pb — 7r(i3) and this equation 
contains the single unknown quantity Pb- Thus, given any 
packet arrival rate A, by solving the equation we can obtain 
Pb, and proceed to derive p-sW by Q. 

Since the relay buffer size is strictly bounded by B (pack¬ 
ets), then no matter what the input rate A is, the relay queue 
is always stable. Eor A G {A : A < ps{X)}, the average delay 
in the local queue is given by flTI 

thus the network is stable. When A ^ {A : A < ps(A)}, for 
the Bernoulli/Bernoulli queue, the queue length will tend to 
infinity. Thus, the network cannot support the input rate stably. 
Since Tc = max{X : X < ps(A)}, then Tc is the throughput 
capacity of the network. ■ 

IV. Simulation 

In this section, we first provide the simulation results to 
validate our theoretical framework, and then we proceed to 
explore how the throughput capacity be influenced by the 
network parameters. 

A. Simulation Setting 

Eor model validation, a simulator in C-H- was developed to 
simulate the packet delivery process in the concerned MANET. 
We focus on a specific node and count its received packets over 
a period of 2 x 10® time slots, to calculate the time-averaged 
throughput. Besides the i.i.d mobility model, another realistic 
mobility model, the random walk model was also implemented 
in the simulator. With random walk mobility model, at the 
beginning of each time slot, each node selects a cell among 
its current cell and its 8 adjacent cells with equal probability, 
then stays in it during this time slot. 

B. Model Validation 

We fix the node number TV = 72, the cell number C = 36, 
and consider three cases of B = 5, i? = 8 and B = 10. 
We increase the packet generation rate A step by step, and 
calculate the service rate ps based on the analytical framework 
provided. The corresponding results are summarized in Eig. 0] 
We can see that for all the three cases there, as A increasing, 
the service rate ps monotonically decreases (note that in the 
infinite buffer scenario, the service rate of local queue has no 
relationship with A ||2l, ifTSl l. and the curves ps = PsiX) and 
















Fig. 4. Service rate fis of the local queue vs. packet generation rate A. 



Fig. 5. Per node’s throughput performance. 


A = A intersect with each other. The intersection points are 
the exact values of throughput capacity for B — 5, B = 8 
and B = 10, respectively. The results are Tc\b =5 = 0.0232, 
Tc\b=8 = 0.0283 and Tc\b=w = 0.0315, it indicates that the 
larger relay buffer size leads to a higher throughput capacity. 

To validate our theoretical results, we summarize the sim¬ 
ulated throughput performance of network scenario {N = 
72, C = 36, i? = 5) in Fig. |5] where the throughput capacity 
Tc is obtained by Fig. |4] and the system load p is dehned as 
the ratio between A and Tc- We can see that the simulated 
throughput linearly increases with p until p = 1, when 
p > 1, the throughput no longer grows up and stays as a 
constant which is consistent with our theoretical analysis. It 
is interesting to observe from Fig. |5] that, the performance 
behavior under random walk mobility model is very similar 
to that under i.i.d mobility model. As shown in El, for a 
cell-partitioned MANET, the throughput capacity under i.i.d 
model is identical to those under non-i.i.d models, if these 
models follow the same-steady distribution. 

C. Performance Analysis 

Fig. |6] shows the relationship between the throughput capac¬ 
ity Tc and the relay buffer size B, under different settings of 



Fig. 6. Throughput capacity Tc vs. relay buffer size B. 



Fig. 7. Throughput capacity Tc vs. number of nodes N. 


node number. The ratio between N and C is hxed by 2. From 
Fig.lSl we can see that for all cases, Tc monotonically increases 
with B, which indicates that the MANET really requires 
a sufficient relay buffer size to guarantee the throughput 
performance. A further careful observation is that when the 
buffer size is large, the increasing of the throughput capacity 
becomes smooth. It indicates that when the buffer size is large 
enough, if we continue to increase it, the performance gain will 
be little. It provides a guideline for the practical applications 
that it is effective to choose an appropriate buffer size to ensure 
the network performance as well as save the cost. 

We proceed to explore the relationship between Tc and 

N, where N/C is hxed by 2. The corresponding results are 
represented in Eig. [T] We can see that for both of the two 
cases there, Tc monotonically decreases with N increasing, 
and vanishes to 0 as tends to inhnity. It is notable that the 
results are quite different from the throughput capacity under 
inhnite buffer scenario, where Tc can keep constant (about 

O. 14 packets/slot) with the network size increasing im, El- It 
indicates that the throughput capacity can not be increased by 
utilizing the node mobility alone, it also requires that the relay 
buffer of each node increases with the network size. 














V. Conclusions 


In this paper, we focus on the throughput capacity of 
MANETs under finite buffer scenario. A modified two-hop re¬ 
lay routing algorithm has been proposed. We have provided an 
ergodic Markov chain-based framework to fully characterize 
the queuing process of relay nodes, and further derived the 
exact throughput capacity. Extensive simulation results have 
been conducted to verify the efficiency of the new theoretical 
framework and show the relationship between Tc and B, N. 
The results indicate that the throughput capacity can not stay 
constant under the finite relay buffer scenario and increasing 
the relay buffer size can upgrade the throughput capacity. 
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